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We discuss the status and the problems related to the application of the off-shell non-perturbative renormal- 
ization method in a fixed gauge to operators relevant to weak decays. In particular, we critically reappraise the 
method recently proposed for the A J =1/2 rule. We also present a general analysis of the renormalization for 
the AJ = 3/2 operators, and apply it to the AS = 2 operator. 



1. Introduction 

Lattice QCD is a unique, systematically im- 
provable, method for computing matrix elements 
from first principles, and has proven a powerful 
and appealing approach. In spite of the successes, 
progress has been slow due to the presence of sys- 
tematic effects, such as discretization and higher- 
order renormalization effects. The Symanzik im- 
provement program is an attrative method 
which allows to reduce discretization order by or- 
der in a in physical quantities. The improvement 
coefficients were computed at first in perturba- 
tion theory (PT) at lowest order, reducing the 
error from 0{a) to ©(^qo) jH^, and recently non- 
perturbatively, achieving a full 0{a?) improve- 
ment In parallel, there has been significant 
progess in the development of non-perturbative 
(NP) methods It is by now generally ac- 

cepted that NP methods in the renormalization 
of lattice operators yield reliable and accurate re- 
sults and should be used whenever possible. In 
the following, we concentrate on the applications 
of the non-perturbative renormalization method 
of ref. 1^, hereafter refered to as NPM. The other 
methods are discussed by Sommer Q and Kura- 
mashi in this Workshop. 

Renormalization of lattice operators is a crucial 
ingredient in the calculation of physical weak ma- 
trix elements on the lattice. A physical amplitude 
Aq_>/3 of a weak transition a — s- /3 is calculated 
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via the Operator Product Expansion (OPE) by 

Ac^^p = Cw{Mw/^l){a\^{^im (1) 

where Cw is the Wilson coefficient of the OPE, 
Mw is the mass of the W boson, ji is the renor- 
malization scale and {a\0{^)\(3) is the matrix el- 
ement of the renormalized operator (at the scale 
^) relevant to the physical process. The Wil- 
son coefficient Cwi^w I contains the short- 
distance information and can be calculated in 
PT in the continuum at the renormalization 
scale /.t. The matrix element contains the long- 
distance dynamics and thus must be calculated 
non-perturbatively on the lattice. Renormaliza- 
tion relates the regularized lattice matrix ele- 
ments to its continuum counterpart. 

On the lattice, chiral symmetry is explicitly 
broken with Wilson-like fermions. The possibility 
of recovering the chiral symmetry in the contin- 
uum limit was shown in ||l3|| . The general pre- 
scription is to subtract from the bare operator 
0{a) all the operators of dimension less or equal 
than 0(a), which have the same quantum num- 
bers conserved by the regularization, 

\im{a\Zo{^xa)[0{a) + y Z,0,{a)m, (2) 

a— >0 ^ — ^ 

i 

where the subtracted operators Oi are not con- 
strained to the same chiral representation of 0{a). 
If these operators have lower dimension than O, 
the mixing constants are power-divergent in the 
cutoff, Z, - l/a'^ with d = dim [O] -dim [O,] > 0. 
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These divergent factors can pick up exponentially 
small contributions in the coupling ag , yielding a 
finite contribution as a — > 0, i.e. 

ie-i/"=W ^Aqcd. (3) 

These divergences must be subtracted in a com- 
pletely NP way. 

2. Operators and phenomenology 

There are a number of four-fermion operators 
which are relevant to different physical processes 
of phenomenological interest. They all have in 
common that their rcnormalization suffers from 
the mixing with chiral violating form factors in- 
duced by the Wilson term. We can divide them 
in three broad classes: 

• LL operators: the A/ = 3/2 components 
are necessary for the Bk parameter, which 
enters in the study of CP violation in K'~' — 
—K^ mixing. It is important to obtain a 
precise and reliable result, as with the mea- 
sured value of the top quark mass, it en- 
ables us to limit the range of values of the 
CP-violation phase S. 

The A/ = 1/2 components on the other 
hand are relevant to the study of the oc- 
tect enhancement in if — > tttt decay, whose 
quantitative understanding still defies the- 
orists. The difference with the A/ = 3/2 
case is that the bare operator are allowed 
to mix with operators of lower dimension- 
ality and hence with coefficients that are 
power-divergent [^jl2|. 

• LR operators: the / — 3/2 part of the 
LR operators which appear in the effective 
weak Hamiltonian due to the electromag- 
netic penguin diagrams are the only oper- 
ators which give rise to an imaginary part 
in the — > tt^tt'^ amplitude, thus yield- 
ing the dominant contribution to e' /e. This 
contribution is usually expressed by the Bj 
and Bs parameters. 

The 1=1/2 part of these operators have 
penguin contractions which make them too 
hard to handle at present. Note also that 



in presence of the heavy top quark, there is 
no GIM suppression. 

• AB — 2 operators: these are the chi- 
ral partners of the A/ = 3/2, and have 
recently been proposed in the study of 
flavour-changing neutral currents in the su- 
persymmetric extensions of the Standard 
Model (m). 

3. The method 

In the NPM, the rcnormalization conditions 
are applied directly to the Green functions of 
quarks and gluons, in a fixed gauge, with given 
off-shell external states of large virtualities 1^. 
The method mimicks what is usually done in the 
perturbative calculation, but the Green functions 
are evaluated in a NP fashion from Monte Carlo 
simulations. 

To give the flavour of the method, let us con- 
sider the simplified case of a multiplicatively 
renormalizable operator, e.g. a two-quark oper- 
ator O — qTq. Given the bare lattice operator 
O^^^^(a), the rcnormalization condition wc im- 
pose is [|j 

= {p\0^^'Ha)\p)\ll%^, (4) 

where {p\ ■ ■ ■ \p) denotes the matrix element of ex- 
ternal quarks of momenta p which can be calcu- 
lated non-perturbatively in the QCD coupling via 
Monte Carlo simulations 1^. The renormalized 
operator obtained with the NPM is then 

6ri(a*) = 4f (5) 

which depends on the external states and the 
gauge, but not on method used to regulate the 
ultra-violet divergences. To stress this point, we 
call the NP rcnormalization scheme Regulariza- 
tion Independent (RI) |l^. The physical opera- 
tor 

OP^y^(M^) = Cni{Mw/ti)Oni{f^) (6) 

is independent of external momenta and gauge 
(up to higher orders in continuum PT and lat- 
tice systematic effects) if the Wilson coefficient 
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function Cuii^^w / in the RI scheme is calcu- 
lated with the same external momenta and gauge 
of ORi(/i). The advantage of the RI scheme is 
that it completely avoids the use of lattice PT, 
which is expected to have a worse convergence 
than the continuum expansion The coeffi- 

cient function CRi(Mvi///i) are instead calculated 
in continuum PT, which cannot be avoided since 
the Wilson OPE is defined perturbatively. 

The phenomenologically more interesting case 
of four-fermion operators are in general not mul- 
tiplicatively renormalizable. The operators which 
need to be subtracted are dictated by the symme- 
tries of the action: charge conjugation (C), parity 
(P) and s <-> c? flavour switching symmmetry (S) 

The main advantage of the NPM is its gen- 
erality, being valid for any composite operator, 
as long as we can can find (a posteriori) a win- 
dow in the range of renormalization scales /x such 
that Aqcd <^ ijl 0(l/a), in order to keep un- 
der control both the higher-order effects in the 
(continuum) perturbative calculation of Crj and 
discretization errors We stress that this re- 
quirement is common to all NP methods on the 
lattice which work at a single value of the lattice 
spacing. The alternative would be matching from 
one value of the coupling to another, as done for 
example in the method of ref . . The main dis- 
advantage of the NPM is the necessity of gauge- 
fixing, which leaves a residual less constraining 
symmetry to dictate the form of the mixing, i.e. 
BRST-invariance. 

4. A/ = 1/2: Seeking new ideas 

A quantitative theoretical understanding of the 
A/ =1/2 rule in if — > tttt decays has proven to be 
a formidable task since the calculation of hadronic 
matrix elements in the low-energy NP regime is 
needed. Let us review the strategies proposed so 
far 0,0. 

In the continuum, with an active charm quark 
and the GIM mechanism at work, the operator 
basis given by 



O 



1 



LL 



-[{sd)L{uu)L ± {su)L{ud)L\ - (m ^ c).(7) 



fermions, the renormalization strategy is com- 
plicated by chiral symmetry breaking. In fact, 
the Wilson term induces the mixing of 
with lower-dimensional operators, with power- 
divergent coefficients, which need to be sub- 
tracted non-perturbatively. In the following, we 
shall concentrate only on the octet component 
of Oj^, which we will denote with . The 
lattice penguin operators, being proportional to 
(m^ — m^a^ <^ 1, will be neglected in the follow- 
ing. 

We are now faced with the problem of 
calculating the four-point matrix elements 
(vr+TT^ |0|^ |iir°). The stardard approach is to rely 
on lowest order chiral PT to relate them to the 
more tractable three-point and two-point matrix 
elements pT| : 

(7r+7r-|0±|i^°) = ^7(±)(™i - m^)//^ (8) 

where ^^^'^ are obtained from 

(7r+(p)|O±lif0(g)) = 7(±)(p . q) - S^^^mllfK 
{Q\0^\K'>) = i5(^\m]^ - ml) 

In this approach, one relies on the calculation 
the K t: matrix element, which only picks up 
a contribution from the parity-conserving (PC) 
part of the operators. Exploiting CPS symme- 
try, we obtain in the PC sector a renormalization 
structure of the form 



O 



PC 



^PC 



Oo± + J2 Zfof + Z^O. + Z^Os 



(9) 



In the framework of lattice QCD with Wilson-likc 



where Oq are the PC bare operators, Of, i — 
1, . . . , 4 are dimension-six operators of wrong chi- 
rality (cf. sec. ||), is the magnetic operator 
sa^^F^^d and O5 is a dimension-three scalar den- 
sity sd. By GIM and power-counting, Z'^ oc 
(m^ — m„) and Z^ oc (m^ — m^jc? . Thus, while 
the coefficient of the magnetic operator can in 
principle be calculated in PT, though it involves 
a two-loop calculation and is very complicated 
fl^ ] , the coefficient of the scalar density is power- 
divergent and can only be reliably calculated in a 
NP fashion. 

There are in principle several NP approaches 
for calculating the mixing coefficients: 
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1. by imposing the Ward Identities on physical 
hadronic states [Tl|] ; 

2. by imposing the Ward Identities on quark 
states, as done for the AS* = 2 operator in 

§; 

3. by the NPM ||. 

The principal drawback of the method 1. is that 
it looses predictive power as the number of coef- 
ficients to determine gets large, as in the present 
case. Methods 2. and 3. are equivalent in the re- 
gion of large /i (cf. sec. ||), and both require the 
inclusion of 

• operators that vanish by the equations of 
motion, because the renormalization condi- 
tions are imposed on off-shell Green func- 
tions. 

Some examples are 

s{p + md)d, s{p + mdfd, (10) 

which, by GIM, generate counterterms propor- 
tional to {rric — mu)/a and (m-c — TOu), respec- 
tively. When inserted in correlations, these terms 
give a finite contribution to the subtraction coef- 
ficient Cs of the scalar density sd. In PT, this 
problem is of course also present, but the spuri- 
ous contributions can be eliminated by looking at 
the momentum dependence 

Cs = C^™° + / + g(j/+i/) + Hp'+p'' + bp ■ p'), 

where /, g, h are calculable functions of and 
rud and p, p' are the momenta of the two external 
legs. In principle, this could also be attempted in 
a NP approach, but it highly unlikely to be able 
to achieve the necessary accuracy. 

• operators are not gauge invariant, because 
the off-shell Green functions are calculated 
in a fixed gauge. 

The non-gauge invariant operators that may mix 
are dictated by lattice BRST symmetry [l9| . 
Some of the these operators have been classified 
in ref. @, e.g. 

smP + md)d, s{P + mdWM (11) 



Again all of these operators need to be taken into 
account as they give a finite contribution to the 
subtraction coefficient. 

The application of the NPM without these ad- 
ditional operators has been outlined in ref. [T^ . 
Let us recall it to give a flavour of its complexity. 
According to the NPM, the mixing Z's are deter- 
mined by finding a set of projectors on the tree- 
level amputated Green functions (GF), with off- 
shell quark and gluon external states, the choice 
of which depends on the nature of the operators 
at hand. For the A/ =1/2 operators we choose 
the following set of external states: qq, qqg, qqqq, 
with the momenta given below in eq. (|l2|). For 
each choice of external states, i.e. for each differ- 
ent set of GF, we need different type of projec- 
tors. Let us denote with Pg the projector on the 
qq GF of the operator O5, with P^ the projec- 
tor on the qqg GF of the operator Oa , and with 
Pj , j = 1, . . . , 4 the set of mutually orthogonal 
projectors on the operators Oi, z = 1, . . . , 4 
Applying the projectors to the corresponding NP 
GF of the renormalized operators O*, with an 
appropriate choice of the external states, we re- 
quire that the renormalized operators be propor- 
tional to the bare operators, 0^{^) oc 0^(a) (up 
to terms of C(a)), i.e. we impose the following 
renormalization conditions (trace over colour and 
spin is understood in the projection operation): 

Ps(9(p)|0±|g(p)) =0 

P^(g(p-%(A:)|6±|g(p)) =0 (12) 
{q{p)q{p)\d^\q{p)q{p)) = 0, j = 1, . . . ,4 

where p and k denote the momentum of the ex- 
ternal quark and gluon legs. The system of equa- 
tions (^2|) in principle completely determines in 
a NP way the renormalization constants, as we 
have six conditions (non-homogeneous due to the 
matrix elements of , cf. eq. (|^)) in six unknown 
mixing constants, , i = 1, . . . , 4, , . 

Unfortunately, since solving eq. (|lj) involves 
delicate cancellations between large contribu- 
tions, it results in a very noisy determination, 
even with large statistics. The main computa- 
tional difficulty lies in the calculation of the GF 
with penguin contractions. The need to include 
the operators which vanish by the equations of 
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motions and are not gauge invariant renders the 
application of the NPM, which was already com- 
plicated without them, highly impractical. We 
conclude that the stardard positive parity meth- 
ods are not a viable way of approaching the 
A/ =1/2 rule and new ideas are needed ||l7| , pT|| . 

5. A/ = 3/2: A general analysis of 
dimension-six mixing 

We now turn to discuss the renormalization 
of A/ = 3/2 operators. They differ from the 
A/ — 1/2 operators in that there are no lower- 
dimensional operators with the same flavour con- 
tent with which they can mix. This implies that 
we need not take into consideration the operators 
that vanish by the equations of motion or the non- 
gauge invariant ones as they can only affect the 
mixing with the lower dimensional operators. 

In order to address this problem, it is conve- 
nient to work with 4 distinct fermion flavours 
V"/, / = 1,...,4, of degenerate mass. Once 
the mixing of the dimension-six generic operators 
with others of the same dimension has been ob- 
tained with four distinct flavours, it is straight- 
forward to apply it to the appropriate operators 
of physical flavours. 

We define the generic four fermion operators 

Or(i)r(^) = (^ir(i)^2)(^3r(2)V'4) 

o,.r(i)t"r(^) = (v;ir(i)rV2)(V53r(2)t°V4) 

of(,,j,<„ = (v;ir(i)V'4)(^3r(2)v.2) 

Of„r(i)t»r(=) = (V^ir(i)iV4)(Vi3r(2)t-^2) 

where F*-^) and F^^^ denotes any Dirac matrix, 
and the colour matrices. Under renormaliza- 
tion, the operators of eq.(p^) can in principle 
mix with any other dimension-six operator, pro- 
vided it has the same quantum numbers. The 
generic QCD Wilson lattice action with 4 de- 
generate quarks is symmetric under parity P, 
and charge conjugation C. Moreover, there are 
three other useful (flavour) symmetries of the 
action, namely the flavour exchange symmetry 
S = (V'2 "04) and the switching symmetries 
S' = (-01 <-> 1/^2, "03 ^ ^^4) and S" = (V'l ^ 
■04, "02 03) In Table |^ we classify the 

operators Op(i)p(2) or combinations of them, ac- 
cording to the discrete symmetries P, C, S' and 



0'r(i)r(2) 
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Ob 


Orb 


Orb 


Ovv 


+1 


+ 1 


+ 1 


+1 


+1 


Oaa 


+1 


+ 1 


+ 1 


+ 1 


+1 


Opp 


+1 


+ 1 


+ 1 


+ 1 


+1 


Oss 


+1 


+ 1 


+ 1 


+ 1 


+1 


Ott 


-t-1 


+ 1 


+ 1 


+ 1 


+1 


0[VA+AV] 


-1 


-1 


-1 


+ 1 


+1 


0[VA-AV] 


-1 


-1 


+ 1 


+ 1 


-1 


0[SP+PS] 


-1 


+ 1 


+ 1 


-1 


-1 


0[sp-ps] 


-1 


+ 1 


-1 


-1 


+1 




-1 


+ 1 


+ 1 


-1 


-1 



Table 1 



Classification of four-fermion operators accord- 
ing to lattice symmetries. These propeties are 
also vahd for the operators Ofay-Wf^T^'^) ■ For the 
operators Op(i)p(2) and 0^p(i)jap(2) : we must ex- 
change the entries of the columns CS' <-> CS" and 
CPS' ^ CPS". 



S". We adopt the notation 

0[r(i)r(2)±r(2)r(i)] = Or(i)r(2) ± Or(2)r<i) 

Note that the results of Table |l] apply also to 
the operators Ojop(i)jap(2) since, upon perform- 
ing the symmetry transformations, sign differ- 
ences, resulting from the presence of the colour 

matrix, disappear because the colour matri- 
ces appear quadratically. On the other hand, 
Op(i)p(2) is obtained by applying S on Op(i)p(2). 
Since S transforms S' into S", the properties of 
Table |l| also apply to Op(i)p(2), but with all S' 
and S columns exchanged. Again, the operator 
O^P(i)jap(2) has the same properties as Op(i)p(2), 
since the colour matrix t"' appears quadratically. 

Our aim is to find complete bases of operators 
which mix under renormalization. Thus, besides 
classifying them according to their symmetries, 
we must also eliminate the operators which are 
not independent. This is seen by applying the 
standard identity of colour matrices 

tAptcD = -^^abScd + ^SadScb (14) 

on the i^'s of a given operator. For the operator 
Ofap(i)jap(i) the result has the general form 

Ojap(,)tap(j) — -TrTy-Op(i)p(j) -I- - ^ C„mOp(„)p(„) 
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where the sum runs over all the Dirac matrices 
obtained by the Fierz transformation of F^^^F^-'^ 
and the factors Cnm are the appropriate con- 
stants of the Fierz transformation |2^. Anal- 
ogously we can express 0^p(i)j„p(j) in terms of 
Or(i)r(j) ^'iid Op(i)p(j). Therefore, in the follow- 
ing, it is adequate to limit ourselves to the mixing 
of Or(»)r(3)'s and Op(i)p(j)'s, according to the en- 
tries of Table 

Having eliminated the non-independent opera- 
tors we proceed in classifying the complete bases 
of operators which mix under renormalization, ac- 
cording to the following two rules: 

1. All operators with identical values of P, 
CPS' and CPS" are allowed to mix with 
each other. These now form a maximal ba- 
sis. 



this basis into 



^[VA+AV] 



0[VA+AV] ± O^^^^y] 



(15) 



and note that Oj^^^^^j has S = +1 and 



^[VA+AV] ^^^^ ^ 



-1. Thus, they do not mix 
with each other. The final result is that the orig- 
inal basis of two operators has been decomposed 
into two bases of one operator each: the two op- 
erators Oj^^^^yj of eq.(p^ renormalize multi- 
plicatively. 

We now turn to O^ya^ay^ . It mixes with 
0^gp_pgY since they both have CPS' -1-1 and 
CPS" = -1. Similarly, Ofy^^^^y^ and 0[sp-ps] 
have CPS' = -1, CPS" = +1. It is convenient to 
combine the two bases into a product basis of 4 
operators: 



2. If possible, the maximal basis must be de- 
composed into smaller bases, by using the 
remaining symmetry S, in order to form lin- 
ear combinations of the operators of our ba- 
sis, which have definite S, i.e. S = -1-1 or 
S = — 1, that can only mix among them- 
selves. 

The first rule is easy to apply in practice, be- 
cause the P, CPS' and CPS" values of the opera- 
tors can be read-off from Table |^ . As an example 
of this rule, we note that 0[gp_pg^ mixes with 
with since they both have P — —1, 

CPS' = -1 and CPS" = -t-1. Having thus ap- 
plied the first rule, we turn to the specific task of 
reducing the basis, for each case of interest (sec- 
ond rule). This we now do case by case, using S 
symmetry. 

5.1. Parity violating operators 

We consider first the parity violating four- 
fermion operators, cf. Tab. |^. None of the four 
violating operators have identical CPS' and CPS" 
values. Each of them, however, mixes with some 
S-counterpart, e.g. 0[sp+ps] and O^p^p^j or 

0[YA-AV], Ofsp_psy 



We examine first O 



O 



[VA+AV] 



which mixes with 



only, forming a basis of two operators, 
characterized by CPS' = CPS" = +1. We rotate 



[VA+AV] 



[VA-AV] 
± 

[SP-PS] ' 



{0[YA-AV] ± 



[SP-PS] 



±0 



[SP-PS] 



(16) 



None of these operators have definite CPS' or 
CPS". However, they have definite S = ±1. 
Thus, they mix in pairs according to their S value; 
i.e. we have reduced the original basis of four op- 
erators into two bases of two operators each. 

Similarly, we rotate the operators 0[sp+ps]j 
0[SP+PS]^ Ott and O^- (with CPS' = CPS" = 
— 1), into the new basis 



^[SP+PS] — 2 



O 



[SP+PS] 



±0 



[SP+PS] 



O: 



TT 



(17) 



which, once more, is decomposed into two bases, 
of two operators each, with definite S = ±1. 

If we introduce, for notational compactness, the 
notation 



Qi — O^YA+AV] 
Q2 = C[yA-^y] 



Qt = 
Q 

the renormalization structure becomes 



3 — ^[SP-PS] 

4 — '-^[SP+PS] 



(18) 



0. 



-5) 



(19) 
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^ 







^23. 













-233 

















Z44 


-2^45 


I 








-254 


-255 / 



where denote the renormaUzed operators and 
Z'^ is the renormahzation (and mixing) matrix 
(summation over repeated indices is impHed). 

Dropping, for simpUcity, the ± subscripts, the 
matrix is a (relatively sparse) block diagonal 
matrix of the form 



(20) 



It is important to notice that this NP renormal- 
ization structure, determined by the symmetries 
of the action, is the same as in the continuum 
nai've dimensional regularization scheme, or any 
other regularization that does not break chirality 
explicitly. In fact, the operators that mix belong 
to the same chiral representation, and their chi- 
ral structures can be obtained from each other by 
Fierz transformations. 

5.2. Parity conserving operators 

Let us now pass to the parity-conserving opera- 
tors, cf. Tab. 1. All of the parity conserving oper- 
ators Orr are eigenstates of all the discrete sym- 
metries listed, with eigenvalue -1-1. Thus by rule 
1., unlike the parity violating case, they all mix 
among each other and also with the five Ofp's; 
the complete maximal basis consists of 10 opera- 
tors. By rule 2., we rotate our basis into a new 
one: 

1 
2 



O 



rr 



[Orr ± O 



^ 1 
rrj 



(21) 



in which the 5 Opp's have being S = -1-1 and 
mix only among themselves; the same is true for 
the Opp's which have S = — 1. Thus the original 
basis of 10 operators has been decomposed into 
two independent bases of 5 operators each. 

This result can be used in the renormahzation 
of the operators 0^y_^^^^, which are the parity 

conserving parts of the operators O^^. Clearly, 

■AA]' ^SS' 



'^PV+AA] ™^eS with 0^y_ 



Orprp. Any other linear combination of these mix- 
ing operators of "wrong" naive chirality is in prin- 
ciple acceptable; particular choices are a question 
of convenience. Here we discuss three such op- 
tions. 



The first option is the one which enables a com- 
parison of this mixing to the perturbative calcu- 
lations of 1 26 1 (and also the NP computations of 
||2^ , |25| , [l^ ). We call this a perturbative-inspired 
(PI) basis: 



Of 



iO± 

4^[VV+AA] 

~'m^'^[ss-pp] 

16Na ^[VV-AA] 



(22) 



o± 

^[SS+PP+TT] 

^o± 

[SS+PP-^TT] 



In this base, Oq mixes with the operators 
of, i — 1, 2, 3 which already appear at the level 
of the one- loop perturbative calculation 1 26 1 , but 
also with of which is not present at the one-loop 
level. The arbitrary numerical overall colour fac- 
tors of of, i — 1, 2, 3 are defined so as to be in 
agreement with the convention of |2^], and the 
colour factor of Of has been set equal to the one 
of of. This choice seems natural for the com- 
parison of its relative weight with respect to the 
other operators present at one-loop. 

A second option, which is exploited in ||] , con- 
sists in taking the basis of the eigenvector of the 
Fierz matrix. We then call this a Fierz-inspired 
(FI) basis: 



Of 



of = o± 



o^ 

^[VV+AA] 



o 



[SS+PP+TT] 

— O^ 

3 — ^[SS+PP-^TT] 

= o^ 

^4 — ^[VV-AA+2{SS-PP)] 
O^ = O^ 

*-^5 — ^[VV-AA-2(SS-PP)] 



(23) 



A third option, used in |^^, consists in the fol- 
lowing basis, which we call Ward-inspired basis, 
as it can obtained, up to signs, from the parity- 
violating base, eq. (jlj), with a chiral transforma- 
tion -04 ^blpi' 



}f 



o 



^ [VV+AA] 
= '^[SS-PP] 



[SS+PP] 
TT 



(24) 



Qf 



±-o± 



o, 
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Whichever basis we choose, upon renormalization 
the structure will be of the form 



,5), 



(25) 



where the matrix is not sparse as in the PV 
case. 



5.3. Scale dependence 

Close to the continuum and chiral limit, the UV 
divergent elements of the renormalization matrix 
Z depend on a/z and only, Zij = Zij{a^, g^), 
whereas non-divergent elements are of the form 
Zij{g^). Since the structure of the renormaliza- 
tion matrix Z is the same as in the continuum, all 
its matrix elements are logarithmically divergent, 
i.e. Zij = Zij{a^i,gl). 

The specification of which elements of the ma- 
trix Z diverge and which are finite, in the limit 
a ~* 0, can be achieved non-perturbatively with 
the aid of the axial Ward Identity (WI). We refer 
to for a detailed presentation while we give 
here only the prescription for the renormalization. 
First, we need to construct the subtracted oper- 
ators 

QT^ = Qi + YTj=2 (^ijQj' 

Qr' = Q. + E,=i4 5CyQj' z = 2,3 (26) 

and, second, renormalize the operators Qf^^ as 
in the continuum, i.e. with a matrix of the form 
given in the PV case, cf. eq. (pO|). 

Before discussing the numerical results, it is 
worth while to stress in which conditions the di- 
rect implementation of the WI, as exploited for 
the AS = 2 operator in (|], and the NPM are 
equivalent. The WI holds for operators with the 
correct chiral properties, that is multiplicatively 
renormalizable operators transforming according 
to a well defined representation of the chiral al- 
gebra. In fact, it is by imposing its validity on 
the renormalized operators one can fix the mixing 
coefficients of the form factors which stem from 
the chiral violation due to the Wilson term. In 
the NPM, this is achieved by imposing that the 
projections of the renormalized operator 0^'S=2 
(cf. eq. (^8|) for the its explicit form) on the four 
chiral violating form factors are zero, 

P,(p|a^^=2|p) = 0, j = l,...,4. (27) 



(a) Z + i 



— 0.2 



(b) Z + 2 



— 0.0 — 
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0.0 



4 - 



— -0.2 — 
5 0.1 



0.5 1 

„2 „2 



Figure 1. NP AS* = 2 mixing Z+'s in PI basis 
as a function of /i^a^ for k = 0.1432. The solid 
(dashed) line is from SPT (BPT). 



But this is true if there are no other causes of 
chiral symmetry breaking, either due to explicit 
presence of mass terms or due to spontaneous 
symmetry breaking in the chiral limit. But both 
these effect die off in the large momenta region. 
So the WI method and the NPM are equivalent 
for sufficiently large values of the renormalization 
scale /i^ = P^- In any case, the overall multiplica- 
tively renormalization constant cancels in the WI 
and thus cannot be determined. Thus the NPM 
(or some other renormalization method) is needed 
even if the WI method is used to obtain the mix- 
ing coefficients. 

6. Numerical results 

Our NP Monte Carlo simulations have been 
performed on an APE machine. We have gen- 
erated an ensemble of 100 independent gauge- 
field configurations, using a tree-level improved 
SW-Clover action action on a 16^ x 32 lat- 
tice, at /3 = 6.0. The quark propagator has 
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(b) z+ 



0.0 
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0.1 
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0,0 — 

— -0.2 — 
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Figure 2. 



dipendence for NP AS* = 2 



mixing Z+'s in PI basis as a function of /i^a^. 
The three symbols correspond to the values of 
K = 0.1425,0.1432,0.1440. 





Q; 




/ 


0.31 


0.030(18) 


0.27(21) 


0.90(15) 


0.62 


-0.027(16) 


0.36(18) 


0.75(13) 


0.96 


-0.012(14) 


0.24(17) 


0.69(12) 


1.27 


0.005(13) 


0.14(16) 


0.68(12) 


1.39 


-0.009(13) 


0.24(16) 


0.67(12) 


1.85 


-0.003(13) 


0.18(16) 


0.66(11) 


2.46 


-0.001(12) 


0.24(15) 


0.65(11) 


4.01 


-0.002(12) 


0.44(15) 


0.67(11) 


BPT 


-0.052(12) 


0.16(15) 


0.62(11) 



Table 2 

Values of a, /3, 7 for 



different values of /ia. 



chosen to present the results at the intermediate 
K = 0.1432. We have chosen not to extrapolate 
in ruqa as we feel that the best one can do with 
a systematic error is to control it rather than ex- 
trapolate in it. 

We first consider the renormalization of the op- 
erator Oq of eq. (p2|), which has the same renor- 
malization properties of parity-conserving part of 
the A5' = 2 operator ajj^dsjl^d. The renormal- 
ized operator is 



been calculated, in the lattice Landau gauge, 
for three values of the hopping parameter k — 
0.1425,0.1432,0.1440, corresponding to a pion 
mass of 900 — 600 MeV. In comparing the NP 
results with PT, we have used both a standard 
bare lattice couphng (SPT) a\^^^ = gl/'^'n: and 
boosted coupling (BPT) of ay = ^^.^"/(Trn) ~ 
1.68 al^" (l|. 

The three values of the hopping parameter al- 
low us to extend the analysis presented in to 
include the study of the mass dependence of the 
Z's, i.e. the effects of the systematics of 0{mqa). 
In general, these effects were expected to be small 
for the light quark sector , though at lower val- 
ues of TOqfl the statistical errors tend to be larger 
and at higher values of mgO the discretization er- 
rors are larger. We have found that indeed the ex- 
pectations are fulfilled, as can be seen from fig. || 
in which the symbols representing the different k 
values can barely be distinguished, at least in the 
significant region of large fi^a^. Thus, we have 



6+ = Z+(/ia)[Oo+(a) + ZtOt{a% 



(28) 



where the operators , . . . , are given 
eq. (E2) for the PI basis, and by 
eq. (p3) for the FI basis. 

In fig. |l| we show the results of the mixing 
Z's in the PI basis at different renormalization 
scales /i^a^. It is clearly notable that that 
and are very well defined and almost scale- 
independent in a large "window" of ^^a^, whereas 
Z^ and Z^ present a smaller window, i.e. a 
more pronounced scale-dependence. Moreover, 
Z^ which is absent in 1-loop PT, is not negli- 
gible. We stress that the large fluctuations at 
small /i^a^ do not spoil the validity of the NPM, 
since in that region the perturbative matching to 
a continuum scheme is not reliable, as for any 
NP lattice method at a fixed lattice spacing 
Fig. ^ shows the mqO-dependence of the Z's in 
the PI basis. As anticipated, the systematic error 
of 0{mqa) is very small, at least in the region of 



m 
in 



10 



Chiral Behaviour 



- <0 + > 






i : 






- i 




'- o SPT □ 


BPT o N = 3 X N = 4 - 







0.05 0,1 0.15 

8/3 f/M/(p.q)/<P5>2 




: 4* 



-o.a 

0.1 



— -0.1 — 



(c) z+. ^ 



-o.a 

0.1 



(b) 



1 0.5 


1 £ 








(d) 











0.5 1 



Figure 3. Chiral behaviour of Bk. 
are taken at js^a^ — 0.96. 



The NP Z's Figure 4. NP A5 = 2 mixing Z+'s in FI basis as 

a function of /i^a^ for k = 0.1432. 



significant vahies of fi^a^, for the quark masses 
we used. 

Using the NP Z+'s in PI basis with their mass 
dependence taken into account (extending the 
analysis of ref. [^), we can revisite the study of 
the chiral behaviour of the Bk parameter [^,^ . 
We adopt the usual parametrization 



(0+) =a + [3m]^ + 7(p • q) + 



(29) 



where a is a lattice artefact that we expect to 
vanish in the chiral limit. In tab. ^ we present 
the parameters a, (3, 7 obtained for different scales 
/i^a^. In fig. H the result using the NP Z+'s at 
/z^ = 0.96 and the bare matrix elements from ||25[| 
is shown. Clearly, the use of the mixing Z's in PT 
(all equal at 1-loop |Q) does not yield the desired 
behaviour, even if a boosted coupling is used. 
This is due to the delicate cancellations which oc- 
cur among the bare matrix elements that can only 
be resolved beyond 1-loop. On the contrary, using 
the NP Z's, the intercept is compatible with zero. 
This behaviour is consistently found at all scale 
/i^a^ > 0.96. Although the use of the complete 



set of operators {N = 4 in fig. yield a better 
chiral behaviour that the use of only the operators 
which mix in 1-loop PT [N = 3), since the mixing 
with O'l starts at O(5o): "^e do not expect dras- 
tic changes in the chiral behaviour. If the chiral 
behaviour were sensibly different, we would not 
trust the matching to the continuum which has 
an uneludable perturbative uncertainty. Indeed 
we find that a^^^ is also compatible with zero, 
and compatible with a^^'* although with large 
statistical errors forced by the thinning approxi- 
mation 1^^. We can only state that the value of 
Bx, proportional to 7, is unaltered and its RGI 
value at NLO is Bk = 0.85 ± 0.15 while the cor- 
rect chiral behaviour of the continuum, signaled 
by the vanishing of a, is recovered |p3| . 

The chiral behaviour of the AS* = 2 has also 
been studied in imposing the chiral WI on 
quark states using the FI basis. The results for 
the mixing Z's, obtained with an unimproved 
Wilson action, seems to show a very stable sig- 
nal as a function of /i^a^ ||l^. In particular, the 
fluctuations at small scales are much reduced. It 
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Figure 5. Overall NP LL Z+ as a function of 
^2a2. The solid (dashed) line is from SPT (BPT). 



Figure 6. Overall NP LL Z as a function of 
ii^a^. The solid (dashed) line is from SPT (BPT). 



could be concluded that the WI method and the 
NPM, although conceptually equivalent at large 
values of /x, are not numerically such. To un- 
derstand this point in more detail, we have tried 
taking linear combinations of the Z's in PI basis 
and expressing them in the FI basis, as shown in 
fig. |4[ It is clear that the Z's, calculated with 
the NPM, in this new basis show a much greater 
stability, comparable to the one obtained by the 
WI method. It must be stressed that although 
the stability of the Z's depends on the choice of 
the basis, the physical results do not. In fact, the 
correct chiral behaviour of the renormalized oper- 
ator is obtained with either basis. This is due to 
the fact that the Z's which present greater fluc- 
tuations multiply bare matrix elements of the op- 
erators that weight less than the ones multiplied 
by the more stable Z's. So that stability of the 
chiral behaviour with the PI basis shown in tab. ^ 
is due to the extremely clean determination of Z2 
in fig. |l| and to the fact that (O2} ~ — 3(Oi 3 4) 

H- 

As a final flourish, in fig. H and fig. we 



show the overall renormalization constants for 

which are the PC and 



Ot,,, , , and 



PV part of the LL operator and, as already 
stressed, can only be determined with the NPM. 
The comparison with PT shows that, while Z+ is 
in good agreement with boosted PT, for Z~ the 
agreement is better with standard PT. 

7. Conclusions 

Recently, there has been considerable progress 
both in the Symanzik improvement program and 
in the development of non-perturbative renormal- 
ization methods. We have presented the applica- 
tion of the off-shell renormalization method using 
a tree-level improved SW-Clover action to four- 
fermion operators with light quarks relevant to 
weak decays. For the A/ = 3/2 sector we have 
presented a general analysis of the renormaliza- 
tion structure, and applied it to the A5 = 2 op- 
erator obtaining the correct chiral behaviour for 
the Bk parameter. On the other hand, for the 
A/ = 1/2 sector, the standard positive parity ap- 
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proach was shown to be unviable, and new ideas 
are needed. 
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